Three-loop corrections to the lightest Higgs scalar boson mass in supersymmetry 



O 

o 

(N 
G 

00 



in 
o 

o 

o 

Oh! 

a: 

CD 
43 



X 



Stephen P. Martin 
Physics Department, Northern Illinois University, DeKalb IL 60115 USA 
and Fermi National Accelerator Laboratory, P0 Box 500, Batavia IL 60510 

I evaluate the largest three-loop corrections to the mass of the lightest Higgs scalar boson in the 
Minimal Supersymmetric Standard Model in a mass-independent renormalization scheme, using ef- 
fective field theory and renormalization group methods. The contributions found here are those that 
depend only on strong and Yukawa interactions and on the leading and next-to-leading logarithms 
of the ratio of a typical superpartner mass scale to the top quark mass. The approximation assumes 
that all superpartners and the other Higgs bosons can be treated as much heavier than the top 
quark, but does not assume their degeneracy. I also discuss the consistent addition of the three-loop 
corrections to a complete two-loop calculation. 
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I. INTRODUCTION 



Low-energy supersymmetry breaking can stabilize the 
electroweak scale against radiative corrections propor- 
tional to much higher mass scales, including the Planck 
mass. In the minimal supersymmetric standard model 
(MSSM), the lightest neutral Higgs scalar boson (h) mass 
is quartically sensitive to the value of the top quark mass, 
but only logarithmically sensitive to the scale of super- 
symmetry breaking, once the Z boson mass is taken as 
fixed. A future experimental determination of the masses 
and couplings of the Higgs scalar bosons and the super- 
partners at the Fermilab Tevatron pp collider, the CERN 
Large Hadron Collider and/or a future linear e + e~ col- 
lider will be crucial in understanding the structure of 
supersymmetry breaking. 

The h mass, in particular, is likely to be a very pre- 
cisely measured quantity [1-4]. This has motivated many 
studies of the relationship between the physical Higgs 
mass Mh and the underlying Lagrangian parameters, in 
the form of radiative corrections of increasing precision 
and detail [5]- [53]. The tremendous effort that has been 
expended on these calculations is necessitated by the ap- 
pearance of qualitatively new enhancement effects at each 
of the first two loop orders in perturbation theory. The 
tree-level result depends only on electroweak gauge cou- 
plings, which enter into the quartic Higgs coupling. At 



one-loop order, the large top Yukawa coupling, enhanced 
by a color factor, enters. At two-loop order, the QCD 
coupling makes its appearance. It turns out that even 
three-loop order contributions will be necessary if the 
goal is to make purely theoretical errors negligible com- 
pared to the future experimental uncertainty in Mh- (Of 
course, there will also be important sources of error due 
to a lack of precise knowledge of the input parameters 
of the theory, such as the top-quark Yukawa couplings 
and the soft supersymmetry breaking terms in the La- 
grangian; these are considered as experimental errors for 
the present discussion.) 

Three general methods have been commonly used, of- 
ten in combination, for evaluating Mh- First, the pole 
mass can be computed by a straightforward calculation 
of the neutral Higgs self-energy diagrams. The resulting 
complete expressions are quite complicated and unwieldy 
beyond one-loop order. A second approach uses the effec- 
tive potential approximation. This means that radiative 
corrections to M% are computed by taking the second 
derivatives of the effective potential; this is equivalent to 
computing the pole mass from self-energy functions in the 
approximation that the external momentum is neglected. 
This has the advantage that calculations can be reduced 
to vacuum graphs, which can always be analytically com- 
puted through 2-loop order. However, this method is not 
gauge-fixing invariant, has limited accuracy, and can suf- 
fer from numerical instabilities if one chooses a renormal- 
ization scale at which a tree-level squared mass happens 
to be extremely small. A third method uses the method 
of effective Lagrangians, with renormalization group run- 
ning used to systematically isolate the effects that are en- 
hanced by logarithms of ratios of the superpartner mass 
scale to the electroweak and top-quark mass scales. Re- 
cent reviews and descriptions of computer programs im- 
plementing some of the known results can be found in 
[29, 46, 49, 50, 53]. 

In this paper, I will use a combination of the three 
methods mentioned above to evaluate the most impor- 
tant 3- loop contributions to Mh- The input parameters 
in this result will be the running DR [54, 55] parame- 
ters in the full theory with no superpartners decoupled. 
(These results can also be converted into on-shell or hy- 
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brid schemes in which all or some of the input particle 
masses are taken to be physical masses rather than run- 
ning masses, although that is not done explicitly here.) 
The results can be used to supplement my previous eval- 
uation of Mh at two-loop order, which includes the full 
diagrammatic results that involve the strong interactions 
and the Yukawa interactions (including the ones that also 
involve the electroweak couplings) [48] , together with all 
other two-loop contributions in the effective potential ap- 
proximation [41, 42]. 

The 3-loop contributions to be found here are only the 
ones that are proportional to powers of the strong cou- 
pling and the top-quark Yukawa coupling. Also, only the 
contributions containing the leading and next-to-leading 
powers of ln(Qg USY /mi) at 3-loop order are evaluated, 
where m t is the top-quark mass and Qsusy is a renormal- 
ization scale comparable to a typical superpartner mass 
scale. However, I will not assume that the superpartners 
are degenerate or that top-squark mixing is negligible. 
These logarithmically enhanced contributions are likely 
to dominate, numerically While there is no unsurmount- 
able obstacle to evaluating the analogous contributions at 
arbitrary loop order using the same methods, as a practi- 
cal matter they are unlikely to be as large as the remain- 



ing uncalculated 3-loop corrections, nor are they likely 
to be as large as the practical experimental errors once 
input parameter uncertainties are taken into account. 



II. CONVENTIONS AND SETUP 



In the following, 



k = 1/16tt 2 



(2.1) 



is used as a loop factor, and the renormalization scale is 
denoted Q. Also, I define the symbol 



ln(ar) = \n(x/Q 2 ). 



(2.2) 



Some formulas below will make use of the one-loop vac- 
uum integral function 



A(x) = x\n(x) — x, 



(2.3) 



and the Passarino-Veltman one-loop self-energy scalar 
loop integral with external momentum invariant s and 
equal internal squared masses x, 



J 



f B (s,x,q 2 ) = 



2 - ln(x/q 2 ) - 2(Ax/s - l) 1 / 2 smT 1 {^/s/Ax) 

2 - ln(x/q 2 ) + (1 - 4.T/s) 1 /2{l n ( s [i _ (1 _ 4x/s) 1 / 2 ]/2x - 1) + m) 

I 



(s < Ax) 
(s > Ax). 



(2.4) 



(This function will appear in some formulas below with 
q not equal to the renormalization scale Q ) This has the 
small s expansion, valid for s < Ax: 

f B (s,x,q 2 ) = -ln(x/q 2 ) + s/(6x) + s 2 /(60x 2 ) 

+s 3 /(A20x 3 ) + ..., (2.5) 

and the special value: 

f B {x,x,q 2 ) = 2-n/VS-Hx/q 2 ). (2.6) 

The neutral Higgs complex scalar field <f> of the Stan- 
dard Model effective theory has a tree-level potential 



V 



X 



H 4 - 



where rh$ and A are running MS parameters. The top 
quark and gauge interactions of the effective Standard 
Model theory are governed by MS running parameters: 



Vt, 33, 9, 9 , v, 



(2.8) 



with other Yukawa couplings neglected. The minimum 
of the tree-level potential occurs at (<f)) — 2fh 2 /X. How- 
ever, here I expand instead around the vacuum expecta- 
tion value (VEV) v defined as the minimum of the loop- 
corrected Landau gauge effective potential of the theory. 



Explicitly, 

Xv 2 = 2ml~i^-AV cS {4>,4> 

v v 0(f> 

For example, working at one-loop order, 



(2.9) 



Xv 2 



2m 2 + K{l2y 2 A(m 2 ) - ^A(m 2 ) 
-i(? 2 +? ,2 )[3A(m|) + 2m|] 
-g 2 [3A(m 2 w ) + 2m 2 w }} +<D(k 2 ), (2.10) 



(2.7) where 



m 2 h 


= Xv 2 , 


(2.11) 


fh t 


= ytv 


(2.12) 


fh.w 


= 9 2 v 2 /2 


(2.13) 


fh 2 z 


- (g 2 + g' 2 )v 2 /2. 


(2.14) 



Equation (2.9) is used to eliminate mi in favor of v. The 
normalization of the Higgs VEV is such that v is roughly 
175 GeV. 

The MSSM theory is governed by (unhatted) running 
DR parameters including the gauge couplings, top-quark 
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Yukawa coupling, and Higgs expectation values (defined 
as the minimum of the loop-corrected Landau gauge ef- 
fective potential of the full MSSM theory): 



53, 9, 9 , Vt, v u , Vd, 



(2.15) 



The last three of these are taken to be real and positive 
by convention. The parameters 



m t = y t v u 
tan(/3) = v u /v d 



(2.16) 
(2.17) 
(2.18) 



arc defined in terms of them, and so depend on the renor- 
malization scale Q. In the following, I will also use the 
short-hand notations 

sp = sin(/3), cp = cos(/3), c 2 p = cos(2/3). (2.19) 

The top-squark sector has a tree-level running squared- 
mass matrix in the (Jl, in) basis: 



ra- 



mi 



■ ■~t v u a t - v dt iy t \ 
v u at - v d ^*Vt m?~ + to 2 J 



(2.20) 



where electroweak £>-terms are neglected (appropriately 
for the approximation used below), and the notation fol- 
lows [41, 56]. The mass eigenstates are related to the 
gauge eigenstates by 



(2.21) 



with |cj| 2 + |sf| 2 = 1. (Here I use the conventions of 
ref. [41]; those of ref. [56] are related by — > — s^.) If yU 
and at are real, then and Cf are the sine and cosine of a 
top-squark mixing angle; otherwise they can be complex 
(but Cj can always be taken real as a convention). It is 
convenient to define the parameter X t by: 

m t X t = v u a t - v d /j,*y t = -s t ~c|(mf 2 - m?), (2.22) 

in terms of which the squared-mass eigenvalues are: 



1 



to? , to? = - [to 2 + to 2 + 2m? 

T{(m 2 R ~ml) 2 +4m 2 \X t \ 2 Y/ 2 ]. (2.23) 



Also appearing below are the tree-level running squared- 
mass eigenvalues of the gluino and squarks, denoted 



2 

mi. 



(* = !,..., 12). 



(2.24) 



The latter include the squared mass eigenvalues of i\ and 
t 2 , as well as the other squarks which are taken to be 
unmixed. The neutralino and chargino mass eigenstates 
will also be taken to be unmixed, with the Higgsinos 
having a common squared mass |/i| 2 . The electroweak 
gauginos do not contribute in the approximation used 
here. Likewise, the Higgs scalar bosons TJ ,^! ,^* are 
treated in the decoupling limit, with a common running 
squared mass m 2 H (supposed to be much larger than m\ 
and to 2 ), and mixing angle a — (i — n/2. 
III. HIGGS POLE MASS IN THE STANDARD 
MODEL 



To prepare for matching the Standard Model to the 
MSSM, one can use the renormalization group to ob- 
tain the higher-loop contributions of leading and next- 
to-leading order in 



1u(Q 2 /to 2 ) 



(3.1) 



for large Q. This is done by using the fact that M% is an 
observable and therefore renormalization-scale indepen- 
dent. Let us write: 



Mi 



(3.2) 



71=0 



p=0 



where the quantities C n . p only depend on Q implicitly 
through the running parameters. For p ^ 0, the coeffi- 
cient C„ iP can be obtained from the results with smaller 
n, provided that the n — p + 1 loop order beta functions 
for each of the running parameters X are known. The 
3-loop order beta functions for the scalar sector of the 
Standard Model are evidently not available at present, 
so only the leading and next-to-leading contributions in 
L can be found at each loop order in this way. In general, 
they satisfy recursion relations: 



d 

dX 



-l,n-l) 



(3.3) 



+ Px ox 



n-2,n-2 



(3.4) 



r 



where 



0x = 



Q^ = ^+n 2 ^ + 



(3.5) 



for X = y t , g 3 , v, A, g, g' . In the following, dependence on 
g, g 1 will be dropped, so the pertinent 2-loop renormaliza- 
tion group equations for the Standard Model parameters 
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are [57, 58]: = -3y 2 v, (3.10) 

& = -24^ + 12Ay 2 + 6A 2 , (3.6) ® = H0? 2 £ 2 + 27^/4 - 3A 2 /8)£, (3.11) 

4? = - 7 »3, (3-12) 



4 2) = -128^ + 120^ + 803 3 2 ^A 

-3^A-36y 2 A 2 -39A 3 /2, (3.7) % = -26g - 2§fy 2 . (3.13) 

4f = 9y 3 /2-8g 2 y t , (3.8) 

^(2) _ _^Qg^4^ t _|_ 3gg 2 y 3 _ 12 j/ 5 ^ n ^ C Standard Model, a routine calculation shows 

5* ' "' ' that the one-loop pole squared mass of the Higgs boson 

— 3yf A + 3y t X 2 /8, (3.9) in the MS scheme can be written as: 

I 



+ \[\{Afh 2 w - ml) - 6g 2 m 2 w }f B (mlrh 2 w ,Q 2 ) + \[\(4m 2 - m\) - 6(g 2 + g' 2 )m 2 z }f B (m 2 h , m 2 , Q 2 ) 



-X[A(m 2 w ) + A(m%)/2} + 2g 2 m 2 w + (g 2 + g 12 )^}. 



(3.14) 



Then, writing C n _ p = c n _ p v 2 7 one can choose: 

co,o - A, (3.15) 
ci, = 2Xy 2 -X 2 (9b h + 2b w + b z + 6/5)/4: 

+0(A 3 ), (3.16) 

from which it follows, via eqs. (3.3)-(3.13), that: 

ci,i = 12^-3y 2 A-3A 2 , (3.17) 
C2,2 = 96g 2 yf-54$-X(l2g 2 y 2 + 99yt/4) 

+ 18A 2 y? + 9A 3 , (3.18) 
C2,i - -32g 2 ^- 18^ + A(20s 2 y 2 + ^[81/20 

-54b h - 12b w - 6b z ]) + 0(X 2 ), (3.19) 
ca,3 = 736^ - 672^ + 90y* + A(-60^ 2 

-lO2g 2 ^ + 243^) + 0(A 2 ), (3.20) 
c 3 ,2 = 160^ + 168^ + [-3246 h - 726^ 

-m z + 6633/10]^ + 0(A), (3.21) 



c M = 5520^- 6492^ + 2178^ 
+783y t lo /2 + 0(A). 



(3.22) 



Here, dependences on g, g' have been dropped except 
where they enter through the kinematic quantities 



2 ~2 ~2\ 
til 



bh = f B (fh h ,m h ,fh 
bw = fB{m 2 h ,fh 2 v ,fn 2 ) 



(3.23) 
(3.24) 

b z = f B (m 2 h ,m 2 z ,m 2 ), (3.25) 
and /s(m 2 ,m 2 ,TO 2 ) has been expanded using eq. (2.5). 



c n , n and c„ i7l _i are needed to generate the terms of order 
X k in the coefficients cn,n and cn.n-i for N > n. How- 
ever, only terms of order A in C n;n and independent of A 
in C n , n -i will be needed in the next section, because A is 
proportional to electroweak couplings at tree level in the 
MSSM. The running parameters in eqs. (3.15)-(3.25) are 
all evaluated at the same arbitrary renormalization scale 
Q appearing in L. 

It should be emphasized that the expansion given 
above for M 2 is far from unique. In particular, the seed 
expression for c^o could have been chosen differently, cor- 
responding e.g. to trading running masses for physical 
masses in the one-loop correction part of eq. (3.14). This 
would produce a different set of higher coefficients, but 
the expression for the total physical mass M% would dif- 
fer only by an amount consistently neglected within the 
approximations. Here, I have chosen to do the expansion 
entirely in terms of running parameters. 



IV. MATCHING TO THE MSSM 

The result for the Higgs pole squared mass M 2 in the 
previous section can now be used to obtain an approxi- 
mate formula in terms of the running DR parameters of 
the MSSM. To do this, one needs the one-loop matching 
conditions, which can be written in the form: 



Note that terms up to order A 



N-n+k 



in the coefficients 



A = 

v = 

Vt = 

93 = 



-(g 2 +g /2 )c% + K yt4c x + ..., (4.1) 

(vn/ap) [1 + ny 2 s 2 c v + ...], (4.2) 

We [1 + k (g 2 3 c vt + y 2 t s 2 c' yt ) + ...], (4.3) 

93 [1 + k.9 2 c 93 +...]. (4.4) 
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The matching coefficients c\, c Vl c Vt , c' yt , and c 93 depend 
on the renormalization scale Q, which is arbitrary but is 
taken to be comparable to the superpartner and heavy 
Higgs bosons (A°, H°, H ± ) masses. All of these scales 
are assumed to be much larger than the top and lightest 
Higgs (h) and electroweak gauge boson masses, so that 
effects suppressed by powers of , m§, mn, etc., 
are neglected. Then one can work consistently to ncxt- 
to-leading order in 

L = ln(Q 2 /m 2 ), (4.5) 



that is, keeping L n and in the terms of n loop order. 
Eliminating the Standard Model parameters in favor of 
MSSM parameters, one obtains: 



n=l 



with 



m 



h 

Ai 
A 2 
A 3 



A 4 



1 l 2 , I2\ 2 2 
2^+5 )C 2 /3V 



12L + c A , 

(96 5 2 - My 2 s 2 )L 2 + [(48c, t - 32). 9 2 + (48^ + 24c„ - 3c A - l8)yls 2 ]L + .... 
(736 ff 4 - 672 ff 2 y 2 S 2 + 90y t 4 4)L 3 + [(160 + 192c S3 + 384c, t ) ff 4 
+ (168 - 12c A - 324c yt + 384^ + I92c v )gly 2 sl 

+ (6633/10- 3246 h - 726^ - 366 z - 99c A /4- 324c^ - 108c v )y^sj]L 2 + . . . , 
(5520g 3 6 - 6492 ff 4 y 2 S 2 + 2178^^4 + 783y^l/2)L 4 + ..., 



(4.7) 

(4.8) 
(4.9) 



(4.10) 
(4.11) 



where dependences on g and g' have been dropped in the loop corrections, except where they enter through the 
kinematic quantities 



b h = f B {ml,ml,rn 2 t ), 
bw = fB{rn\,m^,m 2 ), 
bz = .f B (ml,m 2 z ,m 2 ). 



(4.12) 
(4.13) 
(4.14) 



(For later comparison purposes, the leading- logarithm four loop contribution is also included.) It is important to 
note that the validity of the result just given requires that the expansion is made in terms of running couplings and 
masses, always evaluated at the renormalization scale Q. Indeed, this requirement even applies to the terms that 
are not written here explicitly because they are suppressed by electroweak gauge couplings. [For example, there are 
terms at one- loop order proportional to electroweak couplings multiplied by /^(m^, m|, m 2 ). One could re-express 
those contributions in terms of, for example, /s(M 2 , M§, M 2 ) involving the physical masses M^, Mz and M t , but 
that would require changing the k 3 L 2 coefficient appearing in the expansion above.] 

It remains to find the matching coefficients appearing in the above expressions. First, c\ can be evaluated by 
comparing the well-known one-loop Higgs pole mass calculated directly in the MSSM to eqs. (4.6)-(4.8), giving: 



c\ = 6 



\n(ml) + ln(m? 2 ) + 2\X t \ 2 ln(m? a /m2)/( 



TTlr 
52 

2 /_2 



+ \X t \ 4 {2- [{m 2 2 +m? i )/(m t 2 2 - mf 1 )]ln(<M| 1 )}/K -<) 



(4.15) 



The coefficients c Vt and d + c v are obtained by equating the one-loop expressions for the top-quark pole mass as 
computed in the full MSSM and in the effective Standard Model theory, with the result: 



Cy t + Cy — 



1 + 2/ 1 (m?,m?) + 2/ 1 (m 2 ,m 2 2 ) + 2Rc[X t ]m s {f 2 (m 2 , mfj - / 2 (m?, m?)}/(m? a - m?) 



3 



4(lnm 2 H - 1/2) + fM 2 ,mi) + h(\^,mf 2 ) + /i(|mI^,<) 14 



where 



h(x,y) 



[2a; 2 lna; + 2(y - 2x)y\ny + (x - y)(y - 3x)] /8(x - y) 2 
(ffia;)/4 



{x + 1 y) 
{x = y), 



(4.16) 
(4.17) 

(4.18) 
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f2(x,y) 



1 + y ]n(y/x)/(x - y) [x ^ y) 

{x = y). 



(4.19) 



By comparing the one-loop gluon self-energy functions computed in both the MSSM and the SM effective theory, and 
relying on the equality of physical cross-sections computed in the two theories, one obtains: 



--S3 



(4.20) 



Finally, by comparing the relevant two-loop part of M 2 in eqs. (4.6)-(4.9) above to the known result as calculated 
directly in the MSSM [31, 32, 41, 42], I find: 



-kc t ~| 2 (m? 2 - m\?l(rn\m\) + \X t \ 2 [-3(l - 4| S ^| 2 )m? i m? 2 \^{m\jm\)l1{m\ - m\ 
+3(m| + ml)/ A + \s- t ctf{m\jm\ + m\jm\ - 7m\ - 7mf 2 )/2l /{m\ - m\f . 



(4.21) 



The above results constitute a partial three-loop ap- 
proximation to the lightest Higgs mass in supersymme- 
try. A useful application of this, as has been done earlier 
in [43], is an estimate of the error made in neglecting 
three-loop effects. (See Appendix A for a comparison of 
that paper and others with the results of the present pa- 
per.) However, one would like to go further to use these 
results for an improved calculation of the physical Higgs 
mass. To do so requires consistently adding the three- 
loop correction to a more complete two-loop calculation 
involving electroweak effects, which can be comparable 
in size. 

In earlier work, I have found the two-loop results for 
M 2 in the MSSM, including all diagrammatic contribu- 
tions to the pole mass that involve the strong and Yukawa 
couplings (including those that also involve electroweak 
couplings) [48], as well as all of the remaining contribu- 
tions in the effective potential approach [41, 42]. Since 
the results in the present paper are also given in terms of 
running DR parameters, the three-loop part can be con- 
sistently added to my previous results. However, there is 
an important subtlety involving the identification of the 
tree- level Higgs mass. In refs. [41, 42, 48], the tree- level 
h squared mass is given by the appropriate eigenvalue of 
the (H®,H%) squared mass matrix, evaluated with the 
VEVs at the minimum of the two-loop effective poten- 
tial. In the approximation used here, this corresponds 
to: 



m 2 h 



1 

2^2 



_d_ 

dv v 



+ Vd 



d 
dv d 



AV eS , (4.22) 



where m| was defined by eq. (4.7), and AV e s is the 
radiative part of the effective potential, and the decou- 
pling approximation for the Higgs scalar bosons (rr? h <C 
m 2 H o , m 2 H± , m 2 A o ) has been used. The two versions of the 
tree-level h squared mass, -m 2 h tree and m 2 h , therefore dif- 
fer by tadpole loop contributions that involve the top 
Yukawa coupling and 53. 



To avoid a mismatch between the two-loop part of the 
contribution found in the present paper and the full two- 
loop contribution found in refs. [41, 42, 48], one can take 
the results of those papers and rewrite m\ trec in the tree- 
level and one-loop part in terms of m\ as in eq. (4.22), 
and then expand and incorporate the loop tadpole parts 
as residuals into the one- loop and two-loop parts. In 
the two-loop part, one can consistently simply replace 
m h tree by m\ ■ This is exactly what was done above in 
the derivation of eq. (4.21), albeit in the approximation 
of large gz,yt- (There are no technical obstacles to this 
procedure in general, since the derivatives of the one- 
loop self-energy functions with respect to the external 
momentum invariant and the internal masses are well- 
known, and simple.) The resulting expression, truncated 
at two-loop order, will then allow the three-loop contri- 
bution of eq. (4.10) to be added consistently. 



V. OUTLOOK 

In this paper, I have evaluated the leading and next- 
to-leading logarithm contributions to the lightest Higgs 
mass in the MSSM at three-loop order, in the approxi- 
mation of large QCD and top-quark Yukawa couplings. 
As expected, these contributions are small, but still sig- 
nificant compared to estimates of the experimental er- 
ror for Mh at the LHC or a future linear collider. To 
show the size of the effects, I have plotted the three-loop 
leading- log (L 3 ) contributions in the left panel of figure 
1, in the special limit of a common superpartner mass 
Msusy, and choosing M h = 120 GeV and tan(/3) > 1, 
and using Q = Msusy as the renormalization scale. [See 
eq. (A. 22).] The figure shows the separate contributions 
proportional to a s , as, and independent of as, as well 
as the total. There is a partial cancellation of the a 2 s and 
as parts, which is a fortuitous feature of the perturbative 
expansion scheme chosen here. This illustrates the more 
general fact that the numerical magnitude of the three- 
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1000 1200 



FIG. 1: Leading- logarithm contributions to Mh from 3 loops (left panel) and 4 loops (right panel), as a function of the common 
superpartner mass Msusy- These results are due to the L 3 part of eq. (4.10) and the L 4 part of eq. (4.11), respectively, with 
separate contributions from the various powers of the MSSM QCD coupling and the total shown. The Higgs pole mass is taken 
to be M h = 120 GeV and tan/3 > 1 and af )MS (M z ) = 0.118 (in the five-quark effective Standard Model QCD theory) and 
top-quark pole mass M± = 172 GeV. 



loop correction depends on the way that perturbation 
theory is organized. Changing the renormalization scale 
or scheme, or re-expanding tree-level masses in one-loop 
and two-loop integral functions around pole masses, can 
and does move contributions between loop orders. For 
example, ref. [43] found the result quoted in eq. (A. 17) 
of the present paper, which yields a rather larger numer- 
ical magnitude for the three-loop correction, written in 
terms of Standard Model effective couplings evaluated at 
Q = M t . 

For comparison, the four-loop order leading logarithm 
(L 4 ) contributions are shown in the right panel of figure 
1, using the same scale. As expected, these are quite 
small, and there is again a fortuitous partial cancellation 
between the leading and next-to-leading orders in as- 

The next-to-leading logarithm (L 2 ) contributions at 
three-loop order can be seen to depend on the top-squark 
mixing and other details, and are not depicted numeri- 
cally here, although they can be significant. I expect 
that in application to real-world data, one will want to 
re-express the perturbative expansion by expanding tree- 
level masses in one-loop and two-loop kinematic integrals 
around the pole masses, as this will likely further improve 
the convergence of perturbation theory (see ref. [59] for 
an analogous discussion for the gluino-squark system in 
supersymmctric QCD). This will add more terms to the 
next-to-leading parts of the three-loop contribution. 

If supersymmetry is discovered and thoroughly ex- 
plored at the LHC and a future linear collider, the mass 
of the lightest Higgs boson will present an important pre- 
cision test of our understanding of the theory. The lead- 
ing three- loop corrections will be important for this test. 
It should be emphasized that, at this writing, some two- 
loop contributions remain uncalculated, namely those in- 
volving purely electroweak couplings, which may turn out 



to be similarly important. These contributions cannot be 
captured adequately by the effective potential approxi- 
mation, since some of the relevant self-energy diagrams 
contain a routing by which the external momentum does 
not go through any propagator with mass larger than 
Mz or Mh. Therefore, it will probably be necessary to 
evaluate those two-loop self-energy diagrams on-shell in 
order to reduce theoretical errors to an acceptable level. 

Appendix: Comparison with other work 

It is useful to check the correspondence of the preced- 
ing formulas with earlier special case results. For conve- 
nience, I focus on the two-loop results given in refs. [32] 
and [36] , and the three- loop leading- log results of [43] . 

First, consider the approximations of refs. [32]: 



\s lCi \ = 1/2; (A.l) 

m| i2 = M 2 sT m t \X t \; (A.2) 

ml m t \X t \ « Mf; (A.3) 

m 9 = m i L = m H = H = Ms, (A. 4) 

but \X t \ not necessarily small compared to Mg. Then 
the results above become: 

c vt = l[l-Re[X t ]/M s +E(M|)], (A.5) 

4 = i[(l + c£)E(M|)-4/2], (A.6) 

c x = 12HM 2 s ) + 12\X t \ 2 /M 2 s -\X t \ 4 /Ml, (A.7) 

c v = |At| 2 /4Af s 2 , (A.8) 

c gs = 21n(M|)-l/2, (A.9) 
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in accord with the two-loop results found in eqs. (13)- 
(19) and (25) of ref. [32]. Note that the two-loop non- 
logarithmic parts in eqs. (16) and (17) of ref. [32] are 
not shown explicitly in the present paper; instead, the 
complete contributions up to two- loop order in [41, 42, 
48] can be included by the procedure described at the 
end of section IV. 

Next, consider the approximations of ref. [36], which 
involve a light right-handed top squark: 



mf 



Mf + n 
m t X t 



m t X* 



Ml- 



m 2 , m t \X t \ < M\ < Ml; 



m 5 



rnr 



m H 



M, 



(A.IO) 

(A.11) 
(A.12) 



but \Xt\ not necessarily small compared to Ml. Then 
the results above become 



1 



c Vt = -[l-8Re[Y t ]/M L + 41n(M 



Vt 



c\ = 



6 -(l-^)[l-2HMt)\ 



W 
AMI 



(A.13) 
(A.14) 



61n(A#) + 61n(M^) + — L^- ln(M 2 /M 2 R ) 

' 4 [2 - HMl/Ml )] , (A.15) 

(A.16) 



Mi 



3\X t \ 2 /AM 2 



Li 



in agreement with the two-loop results of Appendix A 
and eqs. (B.8) and (C.12) and (C.13) in ref. [36]. Note 
that ref. [36] also explicitly identifies two-loop contribu- 
tions enhanced by large logarithms ln(M|/M^), which 
are not obtained in the approach used here except when 
they are also enhanced by a ln(M|/raf). This is be- 
cause ref. [36] used a multi-stage effective field theory 
method, first decoupling left-handed top squarks and 
then right-handed top squarks. Also, the approach of 
ref. [36] isolates the non-logarithmic two-loop effective 
potential contributions within the given approximation, 
which are not shown explicitly here. Again, in the ap- 
proach of the present paper, the complete two-loop order 
results of refs. [41, 42, 48] can be included and consis- 
tently supplemented by the three-loop result as described 
at the end of section IV. 



Now, consider the three-loop leading log result given 
in eq. (11) of [43], which assumes a single sparticle mass 
threshold at Ms and tan /3 ^> 1, and is given in terms of 
running MS couplings evaluated at the top mass scale. 
In the notation of section III of the present paper, that 
result is: 



Mi 



where 



12 K L t + K 2 L 2 (-96g 2 + I8y 2 ) 



+^(736^- 240^-99^) 



(A.17) 



L t = ln(M|/m?). (A.18) 
Now, in the leading-logarithm approximation, the trans- 
lation of SM parameters evaluated at the top mass 
scale to the MSSM parameters evaluated at Ms depends 
only the Standard Model one-loop renormalization group 
equations, and is independent of the MS to DR conver- 
sion and threshold corrections. One can write: 



53 = 93 [1 + /sLt(7s|/2) + ...], 
y t = y t [l + KL t (4g 2 -9y 2 /4) 

+n 2 L 2 (22gl - I8g 2 3 y 2 + 243^/32) 



v = 



v[l + K L t (3y 2 /2) 
+ K 2 L 2 (6g 2 y 2 -9yt/4) + ... 



(A.19) 
...],(A.20) 
(A.21) 



where the parameters on the left sides are evaluated at 
Q = m t , and those on the right sides at Q = Ms- Plug- 
ging these into eq. (A.17) immediately yields the special 
form of eqs. (4.6)-(4.10) with only the leading logarithms 
and S/3 = 1: 



Ml = 



4 2 



l2KL t + K 2 L 2 t {96gl~My 2 ) 



+K i Li(736gl- 672giy 2 + 90yt) . (A.22) 



This is more useful as a rough indicator of the sizes of 
theoretical errors due to three-loop effects than as an 
actual precise evaluation of Mh- 

Acknowledgments: This work was supported by 
the National Science Foundation under Grant No. PHY- 
0456635. 



[1] "ATLAS detector and physics performance. Techni- 
cal design report. Vol. 2," CERN-LHCC-99-15, and 
V. Drollinger and A. Sopczak, Eur. Phys. J. C 3, Nl 
(2001), [hep-ph/0102342]. 

[2] T. Abe et al. [American Linear Collider Working Group 
Collaboration] , "Linear collider physics resource book for 
Snowmass 2001. 2: Higgs and supersymmetry studies," 
[hep-ex/0106056]. 

[3] J.A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics 



Working Group Collaboration], "TESLA Technical De- 
sign Report Part III: Physics at an e+e- Linear Collider," 
[hep-ph/0106315]. 
[4] K. Abe et al. [ACFA Linear Collider Working Group Col- 
laboration], "Particle physics experiments at JLC," [hep- 
ph/0109166]. 

[5] SP. Li and M. Sher, Phys. Lett. B 140, 339 (1984). 

[6] M.S. Berger, Phys. Rev. D 41, 225 (1990). 

[7] H.E. Haber and R. Hempfhng, Phys. Rev. Lett. 66, 1815 



9 



(1991). [37 
Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. 
Phys. 85, 1 (1991); Phys. Lett. B 262, 54 (1991). [38 
J. Ellis, G. Ridolfi and F. Zwirner, Phys. Lett. B 257, 83 
(1991); Phys. Lett. B 262, 477 (1991). [39 
R. Barbieri, M. Frigeni and F. Caravaglios, Phys. Lett. 
B 258, 167 (1991). [40 
A. Yamada, Phys. Lett. B 263, 233 (1991). 
J.R. Espinosa and M. Quiros, Phys. Lett. B 266, 389 
(1991). [41 
A. Brignole, Phys. Lett. B 281, 284 (1992). 
M. Drees and M.M. Nojiri, Nucl. Phys. B 369, 54 (1992); 
Phys. Rev. D 45, 2482 (1992). [42 
K. Sasaki, M. Carena and C.E.M. Wagner, Nucl. Phys. 
B 381, 66 (1992). [43 
P.H. Chankowski, S. Pokorski and J. Rosiek, Phys. Lett. 
B 274, 191 (1992); Nucl. Phys. B 423, 437 (1994) [hep- 
ph/9303309]. [44 
H.E. Haber and R. Hempfling, Phys. Rev. D 48, 4280 

(1993) [hep-ph/9307201]. 

J. Kodaira, Y. Yasui and K. Sasaki, Phys. Rev. D 50, [45 
7035 (1994) [hep-ph/9311366]. 

A. Yamada, Z. Phys. C 61, 247 (1994). [46 
R. Hempfling and A.H. Hoang, Phys. Lett. B 331, 99 

(1994) [hep-ph/9401219]. [47 
J. A. Casas, J.R. Espinosa, M. Quiros and A. Riotto, 
Nucl. Phys. B 436, 3 (1995) [Erratum-ibid. B 439, 466 [48 

(1995) ] [hep-ph/9407389]. 

A. Dabelstein, Z. Phys. C 67, 495 (1995) [hep- 
ph/9409375]. 

M. Carena, J.R. Espinosa, M. Quiros and C.E.M. Wag- 
ner, Phys. Lett. B 355, 209 (1995) [hep-ph/9504316]. 
M. Carena, M. Quiros and C. Wagner, Nucl. Phys. B [49 
461, 407 (1996) [hep-ph/9508343]. 

D.M. Pierce, J.A. Bagger, K.T. Matchev and R.J. Zhang, [50 
Nucl. Phys. B 491, 3 (1997) [hep-ph/9606211]. 
H.E. Haber, R. Hempfling and A.H. Hoang, Z. Phys. C [51 
75, 539 (1997) [hep-ph/9609331]. 

S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Rev. D [52 
58, 091701 (1998) [hep-ph/9803277]; Phys. Lett. B 440, 
296 (1998) [hep-ph/9807423]. [53 
R.J. Zhang, Phys. Lett. B 447, 89 (1999) [hep- 
ph/9808299]. [54 
S. Heinemeyer, W. Hollik and G. Weiglein, Com- 
put. Phys. Commun. 124, 76 (2000) [hep-ph/9812320]; 
M. Frank, S. Heinemeyer, W. Hollik and G. Weiglein, [55 
[hep-ph/0202166]; T. Hahn, W. Hollik, S. Heinemeyer 
and G. Weiglein, "Precision Higgs masses with FeynHiggs [56 
2.2," [hep-ph/0507009]. 
[30] S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. J. [57 

C 9, 343 (1999) [hep-ph/98 12472]. 
[31] J.R. Espinosa and R.J. Zhang, JHEP 0003, 026 (2000) 

[hep-ph/9912236]. [58 
[32] J.R. Espinosa and R.J. Zhang, Nucl. Phys. B 586, 3 

(2000) [hep-ph/0003246]. [59 
[33] A. Pilaftsis and C.E.M. Wagner, Nucl. Phys. B 553, 3 

(1999) [hep-ph/9902371]. [60 
[34] M. Carena et al., Nucl. Phys. B 580, 29 (2000) [hep- 
ph/0001002]. 

[35] M. Carena, J.R. Ellis, A. Pilaftsis and C.E.M. Wagner, 

Nucl. Phys. B 586, 92 (2000) [hep-ph/0003180]. 
[36] J.R. Espinosa and I. Navarro, Nucl. Phys. B 615, 82 

(2001) [hep-ph/0104047]. 



Phys. Rev. D 67, 095012 (2003) [hep- 



G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 
611, 403 (2001) [hep-ph/0105096]. 

S. Heinemeyer, Eur. Phys. J. C 22, 521 (2001) [hep- 
ph/0108059]. 

M. Carena, J.R. Ellis, A. Pilaftsis and C.E.M. Wagner, 
Nucl. Phys. B 625, 345 (2002) [hep-ph/0111245]. 
A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. 
Phys. B 631, 195 (2002) [hep-ph/0112177]; Nucl. Phys. 
B 643, 79 (2002) [hep-ph/0206101]. 

S.P Martin, Phys. Rev. D 66, 096001 (2002) [hep- 
ph/0206136]; Phys. Rev. D 65, 116003 (2002) [hep- 
ph/01 11209]. 
S.P. Martin, 
ph/0211366]. 

G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and 
G. Weiglein, Eur. Phys. J. C 28, 133 (2003) [hep- 
ph/0212020]. 

M. Frank, S. Heinemeyer, W. Hollik and G. Weiglein, 
"The Higgs boson masses of the complex MSSM: A com- 
plete one-loop calculation," [hep-ph/0212037]. 

A. Dedes, G. Degrassi and P. Slavich, Nucl. Phys. B 672, 
144 (2003) [hep-ph/0305127]. 

J.S. Lee et al, Comput. Phys. Commun. 156, 283 (2004) 
[hep-ph/0307377]. 

S.P. Martin, Phys. Rev. D 70, 016005 (2004) [hep- 
ph/0312092]. 

S.P. Martin, Phys. Rev. D 71, 016012 (2005) [hep- 
ph/0405022]. The results in this paper are written in 
terms of two-loop integrals reduced to basis integrals 
in S.P. Martin, Phys. Rev. D 70, 016005 (2004) [hep- 
ph/0312092]. The basis integrals are in turn evaluated 
numerically by computer [60]. 

B. C. Allanach, A. Djouadi, J.L. Kneur, W. Porod and 
P. Slavich, JHEP 0409, 044 (2004) [hep-ph/0406166]. 

S. Heinemeyer, Int. J. Mod. Phys. A 21, 2659 (2006) 
[hep-ph/0407244]. 

S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, 
Eur. Phys. J. C 39, 465 (2005) [hep-ph/0411114]. 
W. Hollik and D. Stockinger, Phys. Lett. B 634, 63 
(2006) [hep-ph/0509298]. 

M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak 
and G. Weiglein, [hep-ph/0611326]. 

W. Siegel, Phys. Lett. B 84, 193 (1979); D.M. Capper, 
D.R.T. Jones and P. van Nieuwenhuizen, Nucl. Phys. B 
167, 479 (1980). 

I. Jack et al, Phys. Rev. D 50, 5481 (1994) [hep- 
ph/9407291]. 

S.P. Martin, "A supersymmetry primer," [hep- 
ph/9709356]. 

M.E. Machacek and M.T. Vaughn, Nucl. Phys. B 249, 
70 (1985); Nucl. Phys. B 236, 221 (1984); Nucl. Phys. B 
222, 83 (1983). 

C. Ford, I. Jack and D.R.T. Jones, Nucl. Phys. B 387, 
373 (1992) [hep-ph/0111190]. 

S.P. Martin, Phys. Rev. D 74, 075009 (2006) [hep- 
ph/0608026]. 

S.P. Martin, Phys. Rev. D 68, 075002 (2003) [hep- 
ph/0307101]; S.P. Martin and D.G. Robertson, "TSIL: 
a program for the calculation of two-loop self-energy in- 
tegrals", Comput. Phys. Commun. 174, 133 (2006), [hep- 
ph/0501132]. 



